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INTEGRITY MONITORING FOR GEO-LOCATION SYSTEMS 

FIELD OF THE INVENTION 

The present invention relates generally to wireless communication systems 
and, in particular, to wireless location systems. 

BACKGROUND OF THE RELATED ART 

Satellite-based navigational systems provide accurate, two or three 
dimensional position information to worldwide users. FIG. 1 depicts a well-known 
satellite-based navigational system referred to as Global Positioning System (GPS) 10. GPS 
10 coniprises a plurality of satellites 12-7 and at least one GPS receiver 14. Each sateUite I2-7 
orbits earth at a known speed vj and is a known distance apart from other satellites 12- y. Each 
satellite l2-j transmits a GPS signal 1 1-; which includes a carrier signal with a known 
frequency /modulated using a unique pseudo-random noise (PN-/) code and navigational data 
(ND-f) associated with the particular satellite l2-j, wherein the PN-y code includes a unique 
sequence of PN chips and navigation data ND-j includes a satellite identifier, timing 
information and orbital data, such as elevation angle q, and azimuth angle (j),. 

GPS receiver 14 comprises an antenna 15 for receiving GPS signals 1 a 
plurality of correlators 16-fc for detecting GPS signals 1 l-j and a processor 17 having software 
for determining a position using pseudoranges and the navigation data ND-7 associated with 
the detected GPS signals 1 1-;. GPS receiver 14 detects GPS signals 1 l-j via PN-; codes using 
a correlation process in which correlators 16-fc search for PN-y codes in a carrier frequency 
dimension and a code phase dimension. Such correlation process is implemented as a 
real-time multiplication of a phase shifted replicated PN-y codes modulated onto a replicated 
carrier signal with received GPS signals 11-y, followed by an integration and dump process. 

Upon detecting GPS signals 1 1-y, GPS receiver 14 extracts ranging 
information from the detected GPS signals 1 1-y in the form of pseudoranges or PN phase 
offset measurements, wherein ranging information indicates a range or distance between a 
transmitter source, i.e., satellite 12-y, and a receiver, i.e., GPS receiver 14. The pseudorange 
measurements are subsequentiy used to determine time differences for satellites 12-y 
indicating times required for GPS signals 1 1-y to travel from satellites 12-y to GPS receiver 
14. Such time differences, along with other navigational data indicated in GPS signals 11-y, 
are used to determine the position of GPS receiver 14, as is well known in the art. Note that a 
two dimensional position, i.e., latitude and longitude, of GPS receiver 14 can be determined if 
GPS receiver 14 can detect GPS signals 11-y transmitted from at least three satellites 12-y, and 
a three dimensional position, i.e., latitude, longitude and altitude, of GPS receiver 14 can be 
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determined if GPS receiver 14 can detect GPS signals 1 l-y transmitted from at least four 
satellites ll-j. 

The accuracy of the position being determined for GPS receiver 14 depends 
on the quality of the pseudorange measurements. The quality of the pseudorange 
measurements is affected by signal strengths, multipaths, noise caused by the environment or 
GPS receiver components, etc. Good quality pseudorange measurements result in high 
accuracy GPS receiver 14 positioning solutions. By contrast, poor quality PN phase offset 
measurements, such as outliers, result in GPS receiver 14 positioning solutions with large 
errors. 

Integrity monitoring techniques are used to detect and, if possible, remove 
poor quality pseudorange measurements used in positioning solutions such that high accuracy 
GPS receiver 14 positioning solutions may be obtained. The function of integrity monitoring 
is to perform failure detection and/or failure isolation. Failure detection is a technique for 
determining the existence of a failed satellite, i.e., satellite associated with a poor quality 
pseudorange measurement used in a positioning solution. In order to perform failure 
detection in two dimensional positioning solutions, GPS receiver 14 needs to be able to detect 
at least four GPS signals 1 l-j. In order to perform failure detection in three dimensional 
positioning solutions, GPS receiver 14 needs to be able to detect at least five GPS signals 

11-;. 

Failure isolation is a technique for identifying the failed satellite (or poor 
quality pseudorange measurement). In order to perform failure isolation in two dimensional 
positioning solutions, GPS receiver 14 needs to be able to detect at least five GPS signals 
1 l-j. In order to perform failure isolation in three dimensional positioning solutions, GPS 
receiver 14 needs to be able to detect at least six GPS signals 1 l-j. Once failure isolation is 
successfully performed, the ranging information associated with the failed satellite can be 
removed from positioning solutions. 

Integrity monitoring techniques include the well-known parity method, 
ranging comparison method and least squares residuals method. All three methods are 
snapshot schemes that assume redundant ranging information or measurements, i.e., the 
number of ranging information or measurements is more than required for determining a 
positioning solution, are available at a given sample point in time. All three methods have 
been proven to provide identical results with respect to failure detection. The parity method, 
however, has been proven to be further useful in performing failure isolation. For ease of 
discussion, integrity monitoring will be discussed herein with reference to the parity method. 
This should not be construed to limit the present invention in any manner. 

The basic measurement relationships for failure detection and isolation can be 
generally described by the following well-known equation 
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y = Hx + e equation I 

where}' is a n x 1 measurement vector representing differences between actual measured 
pseudoranges and predicted pseudoranges based on a nominal position of GPS receiver 14 
and a clock bias, n is the number of ranging measurements, and // is a n x m known predictor 
5 matrix arrived at by linearizing about the nominal position and a clock bias. The term jc is a m 
X 1 vector which typically includes components of true position deviation from the nominal 
position plus a derivation of the clock bias, where m is the number of unknown variables in 
vector X to be solved. The term e is a n x 1 measurement error vector due to receiver noise, 
vagaries in propagation, imprecise knowledge of satellite position, satellite clock error, etc. 
10 For the integrity monitoring of a satellite navigation system, it is commonly assumed that e is 
a zero mean with nxn covariance matrix R = o^I, where / is a n x n identity matrix and is a 

^ variance associated with each element of the error vector e . 

in 

.jj In the parity method, a linear transformation on the measurement vector y 

4= results in the following equation 

ru . , ^ . 

P J . 

s where the upper portion yields a least squares solution and the lower portion yields a 

gj parity vector p. Least squares solution jc^ indicates least squares components of true 

position deviation from the nominal position and derivation of the clock bias. Parity vector p 
is the result of operating on measurement vector y with a special (n-m) x n matrix P, wherein 
20 matrix P has rows unity in magnitude and mutually orthogonal to each other and to the 
columns of predictor matrix //. The method of determining matrix P is well-known and 
described in "Matrix Computation," Second Edition, authored by G.H. Golub and C.F. Van 
Loan and published by The John Hopkins University Press. 

Under the assumption of equation 1, i.e., ^ is a zero mean with a n x 1 
25 measurement error vector, failure detection is a simple scalar which obeys a chi-square 

distribution with a freedom of n-m. Specifically, failure detection utilizes a failure decision 
rule, wherein the failure decision rule involves calculating a decision scalar d - p^p^ and 
comparing the decision scalar d against a predetermined failure detection threshold value. If 
the decision scalar d is greater than the failure detection threshold, then failure is detected and 
30 declared. Otherwise, failure is not detected and no failure is declared. Note that n-m must be 
larger or equal to one for failure detection to be possible. Thus, the number of redundant 
ranging measurements n must be at least four or five in order to perform failure detection in 
two or three dimensional positioning solutions, respectively. 

Parity vector p also provides a geometric perspective useful in failure 
35 isolation because the direction indicated by parity vector p can be used to identify the failed 
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satellite. For example, assume there are total 6 satellite measurements and 4 unknown 
variables in vectors, i.e., n=6 and m=4. Assume also a failure happens on a k-th satellite l2-j\ 
which causes an error b in the measurement of the k-ih satellite 12-/ where b is much larger 
than normal measurement noise. The matrix P is then a 2x6 matrix. 



equation 3a 
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Since b is much larger than normal measurement noise, the measurement error vector e can be 
approximated by the following equation. 
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Then, the resultant parity vector p is 



p = Py = P(Hx + €) = Pe== 



Pij 

P2J 



equation 3c 



where p\j and pij denotes elements of matrix P. Equation 3c shows that by ignoring the 
normal measurement noise, parity vector p induced by the error b on satellite 12-; must lie 
along a line whose slope is pi/pir Each satellite 12-y is associated with its own characteristic 
bias line, with a slope determined by the elements of the respective column vector of P. The 
failed satellite is identified as the satellite with a characteristic bias line that lies along parity 
vector p. Note that n-m needs to be greater than or equal to two in order for parity vector p to 
be a vector. If p is not a vector, failure isolation cannot be performed. 

When GPS receiver 14 does not have enough redundant ranging 
measurements available, integrity monitoring can not be performed. For example, GPS 
receiver 14 may not havp a sufficient number of redundant measurements when GPS receiver 
14 is direcdy obstructed from satellites 12-; by a building. In this example, the problem may 
be exacerbated because the number of poor quality pseudorange measurements by GPS 
receiver 14 would probably increase, which would further necessitate the need for integrity 
monitoring. Accordingly, there exists a need to perform integrity monitoring when GPS 
receiver 14 lacks redundant ranging measurements. 



SUMMARY OF THE INVENTION 

The present invention is a method for performing integrity monitoring in die 
absence of redundant ranging measurements from a first ranging source type by using ranging 
measurements from a second ranging source type, wherein a ranging source type corresponds 
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to a class of ranging sources and a ranging source is a transmitter source emitting a signal that 
can be used by a receiver to obtain ranging information indicating a range or distance between 
the transmitter source and receiver. In an embodiment of the present invention, weighted 
ridge regression techniques are used enhance integrity monitoring. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The features, aspects, and advantages of the present invention will become 
better understood with regard to the following description, appended claims, and 
accompanying drawings where 

FIG. 1 depicts a prior art satellite-based navigational system referred to as Global 
Positioning System (GPS); 

FIG. 2 depicts an embodiment of a geo-location architecture comprising of a wireless 
terminal and ranging sources belonging to two different ranging source types in accordance 
with one embodiment of the present invention; 

FIG. 3 depicts a flowchart illustrating a method of integrity monitoring in accordance 
with one embodiment of the present invention; and 

FIG. 4 depicts a flowchart illustrating a method of integrity monitoring in accordance 
with a weighted ridge regression embodiment of the present invention. 

DETAILED DESCRIPTION 

The present invention is a method for performing integrity monitoring of a 
first ranging source type using a second ranging source type. FIG. 2 depicts an embodiment 
of a geo-location architecture 20 comprising of a wireless terminal 22 and ranging sources 
belonging to two different ranging source types, wherein a ranging source type corresponds to 
a class of ranging sources and a ranging source is a transmitter source emitting a signal that 
can be used by wireless terminal 22 (or some other entity, including the ranging source) to 
obtain ranging information indicating a range or distance between the transmitter source and 
wireless terminal 22. In this embodiment, the two ranging source types are a Global 
Positioning System (GPS) 24 and a land based wireless communication network 26. Ranging 
information obtained from wireless communication network 26 will be used to monitor the 
integrity of ranging information obtained from GPS 24, and vice- versa. It should be 
understood that the ranging sources may be a type other than GPS 24 and wireless 
communication network 26. For example, the ranging source type may be another satellite 
system capable of providing ranging information, another wireless communication network 
based on a same or different technology, or some other terrestrial network. It should also be 
understood that a third ranging source may also be used to monitor the integrity of the first 
and/or second ranging source type. 
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GPS 24 comprises a plurality of satellites 28-y, i.e., ranging source. Each 
satellite 28-; orbits earth at a known speed v, and is a known distance apart from other 
satelhtes 28-7. Each sateUite 28-y transmits a GPS signal 30-j which includes a carrier signal 
with a known frequency/modulated using a unique pseudo-random noise (PN-y) code and 
5 navigational data (ND-7) associated with the particular satellite 28-j, wherein the PN-; code 
includes a unique sequence of PN chips and navigation data ND-7 includes a satellite 
identifier, tinning information and orbital data, such as elevation angle a, and azimuth angle 

Ranging information is indicated in GPS signals 3O-7 to wireless terminal 22 
10 in the form of pseudorange or PN phase offset measurements. Wireless terminal 22 includes 
a GPS receiver 32 for measuring pseudoranges or PN phase offsets in detected GPS signals 
3O-7, as is well-known in the art. 

Wireless communication network 26 comprises a plurality of base stations 
34-p, i.e., ranging source. Each base station 34-p transmits a signal 36-p that can be used by 
15 wireless terminal 22 or, with the assistance of wireless terminal 22, base station 34-p to obtain 
nj ranging information (e.g., to obtain round trip delay information). The ranging information 

can be indicated to wireless terminal 22 in a variety of forms. In a first embodiment, ranging 
g information is indicated in the form pilot phase offsets, which correspond to measurements of 

9 code phases in pilot signals. In this embodiment, base stations 34-p are time synchronized 

20 and the pilot signals are transmitted at some known time relative to each other. For example, 
transmission of each pilot signal begins at the same time or at some known time apart. The 
transmitted pilot signals are detected and the pilot phase offsets are measured by wireless 
terminal 22, wherein the pilot phase offset measurements are used to obtain the ranging 
information, as is well-known in the art. 
25 In a second embodiment, ranging information is indicated in the form of 

signal strength. For example, base stations 34-p transmit signals 36-p, such as pilot signals, at 
a known transmission power. Signals 36-p are received and signal strengths measured by 
wireless terminal 22. Based on the known transmission power and the signal strength 
measurements, distances between base stations 34-/? and wireless terminal 22 can be 
30 estimated, as is well-known in the art. 

In a third embodiment, ranging information is indicated in die form of round 
trip delay (RTD) or one way delay (OWD) between base station 34-p and wireless terminal 
22. In this embodiment, base stations 34-p may or may not be time synchronized with respect 
to each other and/or wireless terminal 22. The manners in which RTD and OWD may be 
35 determined are well-known in the art. 



H 3 I 



6 



Da 9-12 



In a fourth embodiment, ranging information is indicated in the form of 
enhanced observed time difference (E-OTD). In E-OTD, unsynchronized base stations 
transmit synchronization and/or dummy bursts which can be used to obtain ranging 
information between the unsynchronized base stations and wireless terminal 22. 

Wireless terminal 22 comprises, along with GPS receiver 32, a correlator 38, 
an oscillator 40 and/or a processor 42 with appropriate software for determining ranging 
information indicated in signal 36-p and for performing integrity monitoring. In one 
embodin^nt, integrity monitoring is performed using one of the following well-known 
techniques: ranging comparison; least square residuals; and parity. All three techniques are 
snapshot schemes that assume redundant ranging measurements, i.e., the number of ranging 
information or measurements is mre than requred for determining a positioning solution, are 
available at a given sample point in time. It has been proven that all three techniques provide 
identical results for failure detection. Among these techniques, the parity technique is 
preferred for failure isolation. For ease of discussion, the present invention will be described 
herein with respect to the parity method. It should not be construed to limit the present 
invention in any manner. 

FIG. 3 depicts a flowchart 300 illustrating a method of integrity monitoring in 
accordance with one embodiment of the present invention. In step 310, wireless terminal 22 
extracts a set of ranging measurements at a given sample point in time from signals 30-; and 
36-p transmitted by satellites 2S-j and base stations 34-/?, respectively. In step 315, the 
extracted ranging measurements are provided as input to processor 42. In step 317, processor 
42 determines if there are redundant ranging measurements sufficient to perform failure 
detection and/or isolation. If no, flowchart 300 return to step 310. If yes, flowchart continues 
to step 320. 

In step 320, processor 42 selects redundant ranging measurements (from the 
set of extracted ranging measurements) to use in the determination of wireless terminal 22's 
position. If two dimensional positioning is being performed by wireless terminal 22, then at 
least five ranging measurements are selected. If three dimensional positioning is being 
performed, then at least six ranging n>easurements are selected. Note that, in step 320, it is 
assumed that both failure detection and failure isolation is being performed as parts of 
integrity monitoring. If only failure detection is being performed, then at least four and five 
ranging measurements are selected for two and three dimensional positioning, respectively. 
Flowchart 300 may include an additional step of determining how many ranging 
measurements to select prior to step 320. For example, it may be determined prior to step 320 
to select only four measurements if there were only four ranging measurements in the set of 
extracted ranging measurements or if only failure detection was to be performed. 
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In a first embodiment, wireless terminal 22 selects all of the extracted ranging 
measurements for determining its position. In a second embodiment, wireless terminal 22 
selects ranging measurements based on the ranging source type. For example, wireless 
terminal 22 selects only ranging measurements extracted from GPS signals 30-7, and not from 
base station signals 36-p, if there were enough ranging measurements extracted from GPS 
signals 30-7 to perform integrity monitoring. In a third embodiment, wireless terminal 22 
selects ranging measurements based on perceived or anticipated reliability of the ranging 
measurements. For example, ranging measurements extracted from base station signals 26-p 
might be selected before ranging measurements extracted from GPS signals 30-) if the signal 
strengths of base station signals 36-p were stronger than the signal strengths of GPS signals 
30-). In a fourth embodiment, wireless terminal 22 selects ranging information from at least 
two different source types. 

In step 330, the selected ranging measurements are inputted into the 
following equation to calculate least squares solution and parity vector p 

\h^H)''H^' 
P 











.P . 





equation 4 



where is a m x 1 vector which typically includes components of true position deviation 
from a nominal position and a clock bias, m is the number of unknown variables in vector 
x^ to be solved, n is the number of ranging measurements, // is a /i x m known predictor 

matrix arrived at by linearizing about a nominal position of wireless terminal 22 and a clock 
bias, and P is a special (n-m) x n matrix with rows unity in magnitude and mutually 
orthogonal to each other and to the columns of predictor matrix H. The method of 
determining matrix P is well-known and described in "Matrix Computation," Second Edition, 
authored by G.H. Golub and C.F. Van Loan and published by The John Hopkins University 
Press. 

The term y represents a n x 1 measurement vector representing differences 
between actual ranging measurements and predicted ranging measurements based on the 
nominal position and clock bias. Measurement vector y is represented by the following 
equation 



y satellite _l2~j 
ybase _ station _ 34- p 



H 



satellite j 



x^- 



^ satellite _n- j 



base _ station _ 34- p 



equation 5 



base ^station JiA-p 

where e is a n x 1 measurement error vector due to receiver noise, vagaries in propagation, 
imprecise knowledge of satellite position, satellite clock error, etc., and the subscripts 
satellite_l2-j and base_station_34'-p indicate the associated satellite 12-y and base station 
34-p for measurement vector y, predictor matrix H and error vector e. 
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In step 340, failure detection is performed by checking if a decision scalar d 
is greater than a failure detection threshold value, where d-p^p. If the decision scalar d is not 
greater than the failure detection threshold value, i.e., failure is not detected, in step 350, the 
least squares solution is outputted as the positioning solution of wireless terminal 22. 
Otherwise, failure is detected and, in step 355, it is determined whether there are enough 
selected ranging measurements to perform failure isolation. If there are enough selected 
ranging measurements, then flowchart 300 proceeds to step 360. Otherwise, flowchart 300 
returns to step 310 via step 365. 

In step 360, failure isolation is performed based on the following equation 



p = Py = Pe=z 



equation 6 



where py and p2j denotes elements of F, e is a n x 1 measurement error vector due to receiver 
noise, vagaries in propagation, imprecise knowledge of satellite position, satellite clock error, 
etc., and yisanxl measurement vector representing differences between actual measured 
pseudoranges and predicted ranges based on a nominal user position and clock bias. The 
failed satellite is isolated as the satellite with a characteristic bias line that lies along parity 
vector p. 

In step 370, ranging measurements are re-selected to be used in determining a 
positioning solution and performing integrity monitoring . In re-selecting the ranging 
measurements, the ranging measurement associated with the failed satellite, if isolated in step 
360, is removed from the group of ranging measurements in which the selection is to be made 
in step 370. In a first embodiment, the re-selected redundant ranging measurements includes 
the remaining ranging measurements, i.e., ranging measurements selected in step 320 less the 
ranging measurement associated with the failed satellite, and previously unselected ranging 
measurements. In a second embodiment, the re-selected ranging measurements includes only 
the remaining ranging measurements if the number of ranging measurements are sufficient to 

r 

perform failure detectfon and/or failure isolation. In a third embodiment, all of the ranging 
measurements are re-selected from a group comprising of the remaining ranging 
measurements and previously unselected ranging measurements in a manner similar to step 
320. Upon re-selecting the redundant ranging measurements, flowchart 300 returns to step 
330. 

Note that the present invention is described herein with integrity monitoring 
being performed at wireless terminal 22. It should be understood that ranging information 
obtained by wireless terminal 22 may also be transmitted to another entity, such as base 
stations 34-p, to have integrity monitoring performed. 
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In another embodiment, the present invention uses weighted ridge regression 
(WRR) techniques to account for different levels of error in ranging measurements related to 
the ranging source type, and for influences of bad geometric distribution. The principles of 
WRR is based on a weighing technique and a ridge regression estimation technique. Errors in 
the ranging information can vary depending on the ranging source type. For example, the 
errors in RTD, pilot phase offset, signal strength and E-OTD measurements are typically 
several times larger than PN phase offset measurements. The weighing technique of WRR 
gives proper weight to each of ranging measurements when determining a positioning 
solution. 

Geometric distribution of wireless terminal 22, base stations 34-p and 
satellites 2S-j differ with various environments. Bad geometric distribution is present when 
the positions of wireless terminal 22, base stations 34-;? and satellites 28-7 located on 
approximately the same plane for three dimensional positioning or the same line for two 
dimensional positioning. In these situations, a position error mechanism known as geometric 
dilution of precision (GDOP) arises. GDOP depends solely on geometry. Position accuracy 
can be estimated as ranging accuracy multiplied by GDOP. Typically, variations in geometry 
are far greater than variations in ranging accuracy. When poor geometry or large GDOP is 
present, position accuracy is adversely affected severely when measurements are applied to 
position solutions by use of the well-known least mean squares (LMS) algorithm. The ridge 
regression estimation technique of WRR minimizes the influence of large GDOP to the mean 
square error of positioning solutions. Weighted ridge regression is described in more detail in 
Appendix A. 

FIG. 4 depicts a flowchart 400 illustrating a method of integrity monitoring in 
accordance with a weighted ridge regression embodiment of the present invention. Steps 410 
to 420 and steps 440 to 470 correspond to steps 310 to 320 and steps 340 to 370, respectively. 
Steps 430 differs from step 330 in that step 430 performs failure detection using the following 
equation 

_P 

where .Jy^ is the weighted ridge regression solution and k is an design parameter in the use 

of the ridge regression method. 

The present invention is described herein with reference to certain 
embodiments. Other embodiments are possible. Accordingly, the present invention should 
not be limited to the embodiments disclosed herein. 
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